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ABSTRACT The engineered control of cellular function through the design of synthetic 
genetic networks is becoming plausible. Here we show how a naturally occurring network 
can be used as a parts list for artificial network design, and how model formulation leads 
to computational and analytical approaches relevant to nonlinear dynamics and statistical 
physics. We first review the relevant work on synthetic gene networks, highlighting the 
important experimental findings with regard to genetic switches and oscillators. We then 
present the derivation of a deterministic model describing the temporal evolution of the 
concentration of protein in a singlc-gene network. Bistability in the steady-state protein 
concentration arises naturally as a consequence of autoregulatory feedback, and we focus on 
the hysteretic properties of the protein concentration as a function of the degradation rate. 
We then formulate the effect of an external noise source which interacts with the protein 
degradation rate. We demonstrate the utility of such a formulation by constructing a protein 
switch, whereby external noise pulses are used to switch the protein concentration between 
two values. Following the lead of earlier work, we show how the addition of a second network 
component can be used to construct a relaxation oscillator, whereby the system is driven 
around the hysteresis loop. We highlight the frequency dependence on the tunable parameter 
values, and discuss design plausibility. We emphasize how the model equations can be used 
to develop design criteria for robust oscillations, and illustrate this point with parameter 
plots illuminating the oscillatory regions for given parameter values. We then turn to the 
utilization of an intrinsic cellular process as a means of controlling the oscillations. We 
consider a network design which exhibits self-sustained oscillations, and discuss the driving 
of the oscillator in the context of synchronization. Then, as a second design, we consider a 
synthetic network with parameter values near, but outside, the oscillatory boundary. In this 
case, we show how resonance can lead to the induction of oscillations and amplification of a 
cellular signal. Finally, we construct a toggle switch from positive regulatory elements, and 
compare the switching properties for this network with those of a network constructed using 
negative regulation. Our results demonstrate the utility of model analysis in the construction 
of synthetic gene regulatory networks. 
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Lead Paragraph 



Many fundamental cellular processes are governed by genetic programs which 
employ protein-DNA interactions in regulating function. Owing to recent tech- 
nological advances, it is now possible to design synthetic gene regulatory net- 
works. While the idea of utilizing synthetic networks in a therapeutic setting is 
still in its infancy, the stage is set for the notion of engineered cellular control at 
the DNA level. Theoretically, the biochemistry of the feedback loops associated 
with protein-DNA interactions often leads to nonlinear equations, and the tools 
of nonlinear analysis become invaluable. Here we utilize a naturally occurring 
genetic network to elucidate the construction and design possibilities for syn- 
thetic gene regulation. Specifically, we show how the genetic circuitry of the 
bacteriophage A can be used to design switching and oscillating networks, and 
how these networks can be coupled to cellular processes. This work suggests 
that a genetic toolbox can be developed using modular design concepts. Such 
advancements could be utilized in engineered approaches to the modification or 
evaluation of cellular processes. 



1 Introduction 

Remarkable progress in genomic research is leading to a complete map of the building blocks 
of biology. Knowledge of this map is, in turn, fueling the study of gene regulation, where 
proteins often regulate their own production or that of other proteins in a complex web of 
interactions. Post-genomic research will likely center on the dissection and analysis of these 
complex dynamical interactions. While the notions of protein-DNA feedback loops and 
network complexity are not new ^ ||, |^, |, |], H, |^, |, ^ |10|, |ll], |12|, experimental advances 



are inducing a resurgence of interest in the quantitative description of gene regulation |13 



14, 15, 16, 17, 18, IS, pO, 21, 22, 23, 24, 25]. These advances are beginning to set the stage for 
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a modular description of the regulatory processes underlying basic cellular function |jT3|, EBl 



27| , |28| , |29| , |34| . In light of nearly three decades of parallel progress in the study of complex 
nonlinear and stochastic processes, the project of quantitatively describing gene regulatory 
networks is timely. 

The concept of engineering genetic networks has roots that date back nearly half a cen- 
tury |31|. It is relatively recent, however, that experimental progress has made the design 
and implementation of genetic networks amenable to quantitative analysis. There are two 
dominant reasons for constructing synthetic networks. First, simple networks represent a 
first step towards logical cellular control, whereby biological processes can be manipulated 
or monitored at the DNA level [^. Such control could have a significant impact on post- 
genomic biotechnology. From the construction of simple switches or oscillators, one can 
imagine the design of genetic code, or software, capable of performing increasingly elabo- 



rate functions [^, |3J] . A second complementary motivation for network construction is the 
scientific notion of reduced complexity; the inherently reductionist approach of decoupling 
a simple network from its native and often complex biological setting can lead to valuable 
information regarding evolutionary design principles . 



Ultimately, we envision the implementation of synthetic networks in therapeutic appli- 
cations. However, such a utilization depends on concurrent progress in efforts to uncover 
basic genomic and interspecies information. For example, broad applicability will only arise 
with detailed information regarding tissue-specific promoters, proteins, and genes. Likewise, 
quantitative network design is contingent on a firm understanding of cellular differentiation 
and fundamental processes such as transcription, translation, and protein metabolism. More 
crucially, delivery is a major hurdle; without identifiable cell-specific recognition molecules, 
there is no method for introducing a network to a specific type of cell. Since, in many re- 
gards, therapeutic applications are somewhat premature, we focus on the implementation of 
synthetic networks in less complicated organisms. The design of synthetic circuits and opti- 
mization of their function in bacteria, yeast, or other plant organisms should reveal nonlinear 
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properties that can be employed as possible mechanisms of cellular control. 

In this paper, we develop several models describing the dynamics of the protein concen- 
tration in small self-contained synthetic networks, and demonstrate techniques for externally 
controlling the dynamics. Although our results are general, as they originate from networks 
designed with common gene regulatory elements, we ground the discussion by considering 
the genetic circuitry of bacteriophage A. Since the range of potentially interesting behavior 
is wide, we focus primarily on the concentration of the A repressor protein. We first show 
how bistability in the steady-state value of the repressor protein can arise from a single-gene 
network. We then show how an external noise source affecting protein degradation can be 
introduced to our model, and how the subsequent Langevin equation is analyzed by way 
of transforming to an equation describing the evolution of a probability function. We then 
obtain the steady-state mean repressor concentration by solving this equation in the long- 
time limit, and discuss its relationship to the magnitude of the external perturbation. This 
leads to a potentially useful application, whereby one utilizes the noise to construct a genetic 
switch. We next show how the addition of a second network component can lead to a genetic 
relaxation oscillator. We study the oscillator model in detail, highlighting the essential de- 
sign criteria. We introduce a mechanism for coupling the oscillator to a time-varying genetic 
process. In the model equations, such coupling leads to a driven oscillator, and we study 
the resulting system in the framework of synchronization. We illustrate the utility of such 
driving through the construction of an amplifier for small periodic signals. Finally, we turn 
to the construction of a genetic toggle switch, and compare switching times for our network 
with those of a network constructed using negative regulation. 



2 Background 



Many processes involving cellular regulation take place at the level of gene transcription . 
The very nature of cellular differentiation and role-specific interaction across cell types im- 
plicates a not yet understood order to cellular processes. Various modeling approaches 



have successfully described certain aspects of gene regulation in specific biological sys- 
tems IP, |12|, |13|, |l^, |18|, ^ ^ ^ It is only recent, however, that designed network 



experiments have arisen in direct support of regulatory models |22| , |23| . In this section, 
we highlight the results of these experimental studies, and set the stage for the discussion of 
the network designs described in this work. 

For completeness, we first discuss the basic concepts of promoters and regulatory feedback 
loops ^ . A promoter region (or, simply, a promoter) denotes a segment of DNA where 
an RNA polymerase molecule will bind and subsequently transcribe a gene into an mRNA 
molecule. Thus, one speaks of a promoter as driving the transcription of a specific gene. 
Transcription begins downstream from the promoter at a particular sequence of DNA that 
is recognized by the polymerase as the start site of transcription. A chemical sequence 
of DNA known as the start codon codes for the region of the gene that is converted into 
amino acids, the protein building blocks. Feedback arises when the translated protein is 
capable of interacting with the promoter that drives its own production or promoters of other 
genes. Such transcriptional regulation is the typical method utilized by cells in controlling 
expression |Q and it can occur in a positive or negative sense. Positive regulation, 
or activation, occurs when a protein increases transcription through biochemical reactions 
that enhance polymerase binding at the promoter region. Negative regulation, or repression, 
involves the blocking of polymerase binding at the promoter region. Proteins commonly 
exist as multi-subunits or multimers which perform regulatory functions throughout the cell 
or serve as DNA-binding proteins. Typically, protein homodimers (or heterodimers) regulate 
transcription, and this fact is responsible for much of the nonlinearity that arises in genetic 
networks |[T9|| . 

Recently, there have been three important experimental studies involving the design of 
synthetic genetic networks. All three employ the use of repressive promoters. In order of 
increasing complexity, they consist of (i) a single autorepressive promoter utilized to demon- 
strate the interplay between negative feedback and internal noise [^, (ii) two repressive 
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promoters used to construct a genetic toggle switch and (iii) three repressive promoters 



employed to exhibit sustained oscillations ||2l|. We now briefly review the key flndings in 
these three studies. 

In the single gene study, both a negatively controlled and an unregulated promoter were 



utilized to study the effect of regulation on variations in cellular protein concentration [23 



The central result is that negative feedback decreases the cell-to-cell fluctuations in protein 
concentration measurements. Although the theoretical notion of network-induced decreased 



variability is not new , this study empirically demonstrates the phenomenon through the 
measurement of protein fluorescence distributions over a population of cells. The flndings 
show that, for a repressive network, the fluorescence distribution is significantly tightened, 
and that such tightening is proportional to the degree to which the promoter is negatively 
controlled. These results suggest that negative feedback is utilized in cellular design as 
a means for mitigating variations in cellular protein concentrations. Since the number of 
proteins per cell is typically small, internal noise is thought to be an important issue, and 
this study speaks to issues regarding the reliability of cellular processes in the presence of 
internal noise. 

The toggle switch involves a network where each of two proteins negatively regulates 
the synthesis of the other; protein "^4" turns off the promoter for gene "i?" , and protein B 
turns off the promoter for gene A In this work, it is shown how certain biochemical 

parameters lead to two stable steady states, with either a high concentration of A (low 
B), or a high concentration of B (low A). Reliable switching between states is induced 
through the transient introduction of either a chemical or thermal stimulus, and shown to be 
significantly sharper than for that of a network designed without co-repression. Additionally, 
the change in fiuorescence distributions during the switching process suggests interesting 
statistical properties regarding internal noise. These results demonstrate that synthetic 
toggle switches can be designed and utilized in a cellular environment. Co-repressive switches 



have long been proposed as a common regulatory theme [ 45 , and the synthetic toggle serves 
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as a model system in which to study such networks. 

In the oscillator study, three repressible promoters were used to construct a network ca- 



pable of producing temporal oscillations in the concentrations of cellular proteins The 
regulatory network was designed with cyclic repressibility; protein A turns off the promoter 
for gene B, protein B turns off the promoter for gene C, and protein C turns off the promoter 
for gene A. For certain biochemical parameters, the "repressilator" was shown to exhibit self- 
sustained oscillations over the entire growth phase of the host E. coli cells. Interestingly, the 
period of the oscillations was shown to be longer than the bacterial septation period, suggest- 
ing that cellular conditions important to the oscillator network were reliably transmitted to 
the progeny cells. However, significant variations in oscillatory phases and amplitudes were 
observed between daughter cells, and internal noise was proposed as a plausible decorrela- 
tion mechanism. These variations suggest that, in order to circumvent the effects of noise, 
naturally-occurring oscillators might need some additional form of control. Indeed, an im- 
portant aspect of this study was its focus on the utilization of synthetic networks as tools 
for biological inference. In this regard, the repressilator work provides potentially valuable 
information pertaining to the design principles of other oscillatory systems, such as circadian 
clocks. 

These studies represent important advances in the engineering-based methodology of 
synthetic network design. In all three, the experimental behavior is consistent with pre- 
dictions which arise from continuum dynamical modeling. Further, theoretical models were 
utilized to determine design criteria, lending support to the notion of an engineering-based 
approach to genetic network design. These criteria included the use of strong constitutive 
promoters, effective transcriptional repression, cooperative protein interactions, and similar 
protein degradation rates. In the immediate future, the construction and analysis of a circuit 
containing an activating control element (i.e., a positive feedback system) appears to be a 
next logical step. 

In this work, we present several models describing the design of synthetic networks in 
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prokaryotic organisms. Specifically, we will utilize genetic components from the virus bacte- 
riophage A. While other quantitative studies have concentrated on the switching properties 



of the A phage circuitry 0, |12| |18|, we focus on its value as a parts list for designing 
synthetic networks. Importantly, the biochemical reactions that constitute the control of A 
phage are very well characterized; the fundamental biochemical reactions are understood, 
and the equilibrium association constants are known ||^, [5^]. In its naturally 

occurring state, A phage infects the bacteria Escherichia coli {E. coli). Upon infection, the 
evolution of A phage proceeds down one of two pathways. The lysis pathway entails the viral 
destruction of the host, creating hundreds of phage progeny in the process. These progeny 
can then infect other bacteria. The lysogenous pathway involves the incorporation of the 
phage DNA into the host genome. In this state, the virus is able to dormantly pass on its 
DNA through the bacterial progeny. The extensive interest in A phage lies in its ability 
to perform a remarkable trick; if an E. coli cell infected with a lysogen is endangered (i.e. 
exposure to UV radiation), the lysogen will quickly switch to the lysis pathway and abandon 
the challenged host cell. 



The biochemistry of the viral "abandon-ship" response is a textbook example ||36| of 
cellular regulation via a naturally-occurring genetic switch. The lytic and lysogenic states 
are controlled by the cro and cl genes, respectively. These genes are regulated by what are 
known as the Prm (c/ gene) and Pr [cro gene) promoters. They overlap in an operator 
region consisting of the three binding sites ORl, 0R2, and 0R3, and the Cro and A repres- 
sor ("repressor", the c/ product) protein actively compete for these binding sites. When 
the Cro protein (product of cro gene) binds to these sites, it induces lysis. When repressor 
binds, lysogeny is maintained and lysis suppressed. When potentially fatal DNA damage is 
sensed by an E. coli host, part of the cellular response is to attempt DNA repair through the 
activation of a protein called RecA. A phage has evolved to utilize RecA as a signal; RecA 
degrades the viral repressor protein and Cro subsequently assumes control of the promoter 
region. Once Cro is in control, lysis ensues and the switch is thrown. 
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3 Bistability in a Single-Gene Network 

In this section, we develop a quantitative model describing the regulation of the Prm operator 
region of A phage. We envision that our system is a DNA plasmid consisting of the promoter 
region and cl gene. 

As noted above, the promoter region contains the three operator sites known as ORl, 
0R2, and 0R3. The basic dynamical properties of this network, along with a categorization 
of the biochemical reactions, are as follows. The gene c/ expresses repressor (Cl), which in 
turn dimerizes and binds to the DNA as a transcription factor. This binding can take place 
at one of the three binding sites ORl, 0R2, or 0R3. The binding affinities are such that, 
typically, binding proceeds sequentially; the dimer first binds to the ORl site, then 0R2, 
and lastly 0R3 Positive feedback arises due to the fact that downstream transcription 
is enhanced by binding at 0R2, while binding at 0R3 represses transcription, effectively 
turning off production and thereby constituting a negative feedback loop. 

The chemical reactions describing the network are naturally divided into two categories 
- fast and slow. The fast reactions have rate constants of order seconds, and are therefore 
assumed to be in equilibrium with respect to the slow reactions, which are described by rates 
of order minutes. If we let X, X2, and D denote the repressor, repressor dimer, and DNA 
promoter site, respectively, then we may write the equilibrium reactions 



x + x 




X2 


D + X2 












Di+X2 




D2D1 


D2D1 + X2 




D3D2D 



where Di denotes dimer binding to the ORz site, and the Ki = ki/k^i are equilibrium 
constants. We let = (J1K2 and = o"2-ft'2, so that ai and (72 represent binding strengths 
relative to the dimer-ORl strength. 

The slow irreversible reactions are transcription and degradation. If no repressor is bound 
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to the operator region, or if a single repressor dimer is bound to ORl, transcription proceeds 
at a normal unenhanced rate. If, however, a repressor dimer is bound to 0R2, the binding 
affinity of RNA polymerase to the promoter region is enhanced, leading to an amplification of 
transcription. Degradation is essentially due to cell growth. We write the reactions governing 
these processes as 

D + p'^D + P + nX (2) 
Di + P Di + P + nX 
D2Di + P D2Di + P + nX 
X h 

where P denotes the concentration of RNA polymerase, n is the number of repressor proteins 
per mRNA transcript, and a > 1 is the degree to which transcription is enhanced by dimer 
occupation of 0R2. 

Defining concentrations as our dynamical variables, x = [X], X2 = [X2], do = [D], 
di = [Di], d2 = [D2D1], and d^ = [D3D2D1], we can write a rate equation describing the 
evolution of the concentration of repressor, 

X = —2kix^ + 2A;_iX2 + nktpoido + di + ac/2) — k^x (3) 

where we assume that the concentration of RNA polymerase Pq remains constant during 
time. 

We next eliminate X2 and the di from Eq. (^) as follows. We utilize the fact that the 
reactions in Eq. (|l|) are fast compared to expression and degradation, and write algebraic 
expressions 

X2 = Kix" (4) 
di = K2doX2 = {KiK2)dox'^ 
d2 = K'idiX2 = ai{KiK2fdQX^ 
ds = K4d2X2 = aia2{KiK2)^dox^ 
10 



Further, the total concentration of DNA promoter sites dx is constant, so that 

mdT = do{l + K^K^x^ + aiiK^K^fx^ + aia2{K^K2fx^) (5) 

where m is the copy number for the plasmid, i.e., the number of plasmids per cell. 

We next eliminate two of the parameters by rescaling the repressor concentration x 



and time. To this end, we define the dimensionless variables x = x^ K1K2 and t 



t{ktPodTny/ K1K2) . Upon substitution into Eq. (^, we obtain 



m(l + + aai x"^) 
X = ^ z T — X (6) 



where = k^/ {dxnktPQyJ Ki K2 ) , the time derivative is with respect to t, and we have 
suppressed the overbar on x. The equilibrium constants are Ki = 5.0 x lO^M^^ and 
K2 = 3.3 X lO^M^^ so that the transformation from the dimensionless vari- 



able X to the total concentration of repressor (monomeric and dimeric forms) is given by 
[CI]= (7.7x + 3.0a;^) nM. The scaling of time involves the parameter kt, and since transcrip- 
tion and translation are actually a complex sequence of reactions, it is difficult to give this 



lump parameter a numerical value. However, in Ref. |5§], it is shown that, by utilizing a 
model for the lysogenous state of the A phage, a consistency argument yields a value for the 
product of parameters {dtnktPo) = 87.6 nM min^^. This leads to a transformation from the 
dimensionless time t to time measured in minutes of t{minutes) = 0.089t. 

Since equations similar to Eq. (j^) often arise in the modeling of genetic circuits (see 
Refs. [Q) of this Focus Issue), it is worth noting the specifics of its functional form. The first 
term on the right hand side of Eq. (^) represents production of repressor due to transcription. 
The even polynomials in x occur due to dimerization and subsequent binding to the promoter 
region. As noted above, the ai prefactors denote the relative affinities for dimer binding to 
ORl versus that of binding to 0R2 (o"i) and 0R3 ((T2)- The prefactor a > 1 on the x^ term 
is present because transcription is enhanced when the two operator sites ORl and 0R2 are 
occupied (x^x^). The x^ term represents the occupation of all three operator sites, and arises 
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in the denominator because dimer occupation of 0R3 inhibits polymerase binding and shuts 
off transcription. 

For the operator region of A phage, we have (Xi ~ 2, 0-2 ~ 0.08, and a ~ 11 P, ^ ^ 



so that the parameters 7^ and m in Eq. (j^) determine the steady-state concentration of 
repressor. The parameter 7^; is directly proportional to the protein degradation rate, and 
in the construction of artificial networks, it can be utilized as a tunable parameter. The 
integer parameter m represents the number of plasmids per cell. While this parameter is not 
accessible during an experiment, it is possible to design a plasmid with a given copy number, 
with typical values in the range of 1-100. 

The nonlinearity of Eq. (^) leads to a bistable regime in the steady state concentration of 
repressor, and in Figure lA we plot the steady-state concentration of repressor as a function 
of the parameter 7^.. The bistability arises as a consequence of the competition between the 
production of x along with dimerization and its degradation. For certain parameter values, 
the initial concentration is irrelevant, but for those that more closely balance production and 
loss, the final concentration is determined by the initial value. 

Before turning to the next section, we make one additional observation regarding the 
synonymous issues of the general applicability of a synthetic network and experimental mea- 
surement. In experimental situations, a Green Fluorescent Protein (GFP) is often employed 
as a measurement tag known as a reporter gene. This is done by inserting the gene encoding 
GFP adjacent to the gene of interest, so that the reporter protein is produced in tandem 
with the protein of interest. In the context of the formulation given above, we can generalize 
Eq. to include the dynamics of the reporter protein, 

X = f{x) - 7x a; (7) 

9 = f{x) -Ig 9 



where f{x) is the nonlinear term in Eq. (^, 7^ = kg/ {dTnktPoV K1K2) , and the GFP con- 



centration is scaled by the same factor as repressor {g = g^/KiK^). In analogy with the 
equation for x, kg is the degradation rate for GFP, and we have assumed that the number 
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of proteins per transcript n is the same for both processes. This abihty to co-transcribe two 
genes from the same promoter and transcribe in tandem has two important consequences. 
First, since proteins are typically very stable, it is often desirable to substantially increase 



their degradation rate in order to access some nonlinear regime Such a high degra- 

dation rate typically will lead to a low protein concentration, and this, in turn, can induce 
detection problems. The utilization of a GFP-type reporter protein can help to mitigate this 
problem, since its degradation rate can be left at a relatively low value. Second, and perhaps 
more importantly, are the significant implications for the generality of designer networks; in 
prokaryotic organisms, any protein can be substituted for GFP and co-transcribed, so that 
one network design can be utilized in a myriad of situations. 

4 A Noise-Based Protein Switch 

We now focus on parameter values leading to bistability, and consider how an external noise 
source can be utilized to alter the production of protein. Physically, we take the dynamical 
variables x and g described above to represent the protein concentrations within a colony of 
cells, and consider the noise to act on many copies of this colony. In the absence of noise, 
each colony will evolve identically to one of the two fixed points, as discussed above. The 
presence of a noise source will at times modify this simple behavior, whereby colony-to-colony 
fluctuations can induce novel behavior. 

Noise in the form of random fluctuations arises in biochemical networks in one of two 
ways. As discussed elsewhere in this Focus Issue , internal noise is inherent in biochem- 



ical reactions, often arising due to the relatively small numbers of reactant molecules. On 
the other hand, external noise originates in the random variation of one or more of the 
externally- set control parameters, such as the rate constants associated with a given set of 
reactions. If the noise source is small, its effect can often be incorporated post hoc into the 
rate equations. In the case of internal noise, this is done in an attempt to recapture the lost 
information embodied in the rate-equation approximation. But in the case of external noise, 
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one often wishes to introduce some new phenomenon where the details of the effect are not 
precisely known. In either case, the governing rate equations are augmented with additive 
or multiplicative stochastic terms. These terms, viewed as a random perturbation to the 
deterministic picture, can induce various effects, most notably, switching between potential 
attractors (i.e., fixed points, limit cycles, chaotic attractors) |]5^ . 

In previous work, the effects of coupling between an external noise source and both 
the basal production rate and the transcriptional enhancement process were examined |[55|| . 
Here, we analyze the effect of a noise source which alters protein degradation. Since the 
mathematical formulation is similar to that of Ref. our goal here is to reproduce the 
phenomenology of that work under different assumptions. As in Ref.|5^, we posit that the 
external noise effect will be small and can be treated as a random perturbation to our existing 
treatment; we envision that events induced will be interactions between the external noise 
source and the protein degradation rate, and that this will translate to a rapidly varying 
protein degradation embodied in the external parameters and 7^. In order to introduce 
this effect, we generalize the model of the previous section such that random fluctuations 
enter Eq. (0) multiplicatively, 

i = fix)-{r,-Ut))x (8) 
9 = f{x)-{7,-Ut))9 (9) 

where the C,i{t) are rapidly fluctuating random terms with zero mean (< ^i{t) >= 0). In order 
to encapsulate the independent random fluctuations, we make the standard requirement 
that the autocorrelation be "(5-correlated" , i.e., the statistics of the C,i{t) are such that < 
C,i(t)C,j(t') >= D6ij{t — t'), with D proportional to the strength of the perturbation, and we 
have assumed that the size of the induced fluctuations is the same for both proteins. 

Since, in Eqs. (|]) and @, the reporter protein concentration g does not couple to the 
equation for the repressor concentration, the qualitative behavior of the set of equations 
may be obtained by analyzing x. We first define a change of variables which transforms the 
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multiplicative Langevin equation to an additive one. Letting x = e^, Eq. (||) becomes, 

1 + + 22e^" 



-7x + ex(t) (10) 



Eq. ([10| ) can be rewritten as: 

50(2) 



^2; 

where the potential (j){z) is introduced: 



+ e.(t) (11) 



<i){z) = - j g{z)dz (12) 

0(2;) can be viewed as an "energy landscape", whereby z is considered the position of a 
particle moving in the landscape. One such landscape is plotted in Fig. IB. Note that the 
stable fixed points correspond to the minima of the potential in Fig. lb, and the effect of 
the additive noise term is to cause random kicks to the particle (system state point) lying 
in one of these minima. On occasion, a sequence of kicks may enable the particle to escape 
a local minimum and reside in a new valley. 

In order to analyze Eq. (pT]), one typically introduces the probability distribution P{z^ t), 
which is effectively the probability of finding the system in a state z at time t. Then, given 



Eq. (|TT]), a Fokker-Planck (FP) equation for P{z, t) can be constructed |5^. The steady-state 
solution for this equation is given by 

Ps{z) = Ae-i^^'^ (13) 

where A is a normalization constant determined by requiring the integral of Ps{z) over all z 
be unity. 

Using the steady-state distribution, the steady-state mean (ssm), < z >ss, is given by 



00 

2 



^^e-u<^W^^ (14) 
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In Fig. IC, we plot the ssm value of z as a function of D, obtained by numerically integrating 
Eq. (II^). It can be seen that the ssm of z increases with D, corresponding to the increasing 
likelihood of populating the upper state in Fig. IB. 

Figure IC indicates that the external noise can be used to control the ssm concentration. 
As a candidate application, consider the following protein switch. Given parameter values 
leading to the landscape of Fig. IB, we begin the switch in the "off" position by tuning the 
noise strength to a very low value. This will cause a high population in the lower state, and 
a correspondingly low value of the concentration. Then at some time later, consider pulsing 
the system by increasing the noise to some large value for a short period of time, followed 
by a decrease back to the original low value. The pulse will cause the upper state to become 
populated, corresponding to a concentration increase and a flipping of the switch to the "on" 
position. As the pulse quickly subsides, the upper state remains populated as the noise is 
not of sufficient strength to drive the system across either barrier (on relevant time scales). 
To return the switch to the off position, the upper-state population needs to be decreased to 
a low value. This can be achieved by applying a second noise pulse of intermediate strength. 
This intermediate value is chosen large enough so as to enhance transitions to the lower 
state, but small enough as to remain prohibitive to upper-state transitions. 

Figure ID depicts the time evolution of the switching process for noise pulses of strengths 
D = 1.0 and D = 0.1. Initially, the concentration begins at a level of ~ 0.4 /iM, correspond- 
ing to a low noise value of D = 0.01. At 40 minutes, a noise burst of strength D = 1.0 is 
used to drive the concentration to a value of ~ 2.2 yuM. Following this burst, the noise is 
returned to its original value. At 80 minutes, a second noise burst of strength D = 0.1 is 
used to return the concentration to its original value. 

5 A Genetic Relaxation Oscillator 

The repressillator represents an impressive step towards the generation of controllable in vivo 
genetic oscillations. However, there were significant cell-to-cell variations, apparently arising 
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from small molecule number fluctuations |21, p5|. In order to circumvent such variability. 



the utilization of hysteresis-based oscillations has recently been proposed |3^. In this work 



it was shown how a model circadian network can oscillate reliably in the presence of internal 
noise. In this section, we describe an implementation of such an oscillator, based on the 
repressor network of Section 3. 

The hysteretic effect in Fig. lA can be employed to induce oscillations, provided we can 
couple the network to a slow subsystem that effectively drives the parameter 73,. This can 
be done by inserting a repressor protease under the control of a separate Prm promoter 
region. The network is depicted in Fig. 2A. On one plasmid, we have the network of Section 
3; the repressor protein CI, which is under the control of the promoter Prm, stimulates its 
own production at low concentrations and shuts off the promoter at high concentrations. 
On a second plasmid, we again utilize the Prm promoter region, but here we insert the 
gene encoding the protein RcsA. The crucial interaction is between RcsA and CI; RcsA 
is a protease for repressor, effectively inactivating its ability to control the Prm promoter 
region 
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The equations governing this network can be deduced from Eq.(|^) by noting the following. 
First, both RcsA and repressor are under the control of the same promoter, so that the 
functional form of the production term f{x) in Eq.(|^) will be the same for both proteins. 
Second, we envision our network as being constructed from two plasmids - one for repressor 
and one for RcsA, and that we have control over the number of plasmids per cell (copy 
number) of each type. Lastly, the interaction of the RcsA and repressor proteins leads to the 
degradation of repressor. Putting these facts together, and letting y denote the concentration 
of RcsA, we have 

X = m^f{x)--f:rx-jxyxy (15) 

= m^fix) --f{y)x 
y = myf{x) - -fyy 
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where 7(7/) = 7^ + jxyV, and rrix and rUy denote the plasmid copy numbers for the two 
species. 

In Fig. 2B, we present simulation results for the concentration of repressor as a function 
time. The nature of the oscillations can be understood using Fig. lA. Suppose we begin 
with a parameter value of 7(1/) = 4 and on the upper branch of the figure. The large 
value of repressor will then serve to activate the promoter for the RcsA, and thus lead to 
its increased production. An increase in the RcsA acts as an additive degradation term 
for repressor (see Eq. |T5|), and thus effectively induces slow motion to the right on the 
upper branch of Fig. lA. This motion will continue until the repressor concentration falls 
off the upper branch at 7(1/) ~ 5.8. At this point, with the repressor concentration at a 
very low value, the promoters are essentially turned off. Then, as RcsA begins to degrade, 
the repressor concentration slowly moves to left along the lower branch of Fig. lA, until it 
encounters the bifurcation point at 7(7/) ~ 3.6. It then jumps to its original high value, with 
the entire process repeating and producing the oscillations in Fig. 2B. 

The oscillations in Fig. 2B are for specific parameter values; of course, not all choices 
of parameters will lead to oscillations. The clarification of the specific parameter values 
leading to oscillations is therefore important in the design of synthetic networks ||21|. For 
proteins in their native state, the degradation rates and 'jy are very small, corresponding 
to the high degree of stability for most proteins. For example, a consistency argument 
applied to a similar model for A phage switching leads to 7^ ~ 0.004. However, using 
a temperature-sensitive variety of the repressor protein, 73. can be made tunable over many 
orders of magnitude. Other techniques, such as SSRA tagging or titration, can be employed 
to increase the degradation rate for RcsA. The copy numbers rrix and rriy can be chosen for a 
particular design, and the parameter •y^y, which measures the rate of repressor degradation 
by RcsA, is unknown. 

In Fig. 3A, we present oscillatory regimes for Eq. (p!5[) as a function of •jx and ■jy, and 
for two fixed values of the parameter We see that the oscillatory regime is larger for 
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smaller values of the parameter ^^y. However, the larger regime corresponds to larger values 
of the degradation rate for RcsA. Interestingly, if we take the native (i.e., without tuning) 
degradation rates to be 7a; ~ 7j/ ~ 0.005, we note that the system is naturally poised very 
near the oscillatory regime. In Fig. 3B, we present the oscillatory regime as a function of 
7x and 73:2;, and for two fixed values of 7^. The regime is increased for smaller values of 
7j^, and, in both cases, small values of 7^; are preferable. Moreover, both Figs. 3A and 3B 
indicate that the system will oscillate for arbitrarily small values of the repressor degradation 
parameter 7^;. In Fig. 3C we depict the oscillatory regime as a function of the copy numbers 
rUx and m^, and for fixed degradation rates. Importantly, one can adjust the periodic regime 
to account for the unknown parameter ^^y- Figure 3C indicates that, for oscillations, one 
should choose as large a copy number as possible for the plasmid containing the repressor 
protein {rrix)- Correspondingly, one should design the RcsA plasmid with a significantly 
smaller copy number rUy. 

We now turn briefly to the period of the oscillations. If designed genetic oscillations are to 
be utilized, an important issue is the dependence of the oscillation period on the parameter 
values. In Fig. 4A we plot the oscillation period for our CI-RcsA network as a function of 
the degradation parameter and for other parameter values corresponding to the lower 
wedge of Fig. 3A. We observe that an increase in will decrease the period of oscillations. 
Further, since the cell-division period for E. coli is ~ 35 — 40 minutes, we note that the lower 
limit roughly corresponds to this period, and that, at the upper limit, we can expect four 
oscillations per cell division. The utilization of tuning the period of the oscillations to the 
cell-division time will be discussed in the next section. In Fig. 4B, we plot the period as a 
function of the copy number rrix- We observe that the period depends very weakly on the 
copy number. 
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6 Driving the Oscillator 

We next turn to the utilization of an intrinsic cellular process as a means of controlling the 
oscillations described in the previous section. We will first consider a network design which 
exhibits self-sustained oscillations (i.e., with parameters that are in one of the oscillatory 
regions of Fig. 3), and discuss the driving of the oscillator in the context of synchronization. 
As a second design, we will consider a synthetic network with parameter values near, but 
outside, the oscillatory boundary. In that case, we will show how resonance can lead to the 
induction of oscillations and amplification of a cellular signal. 

We suppose that an intrinsic cellular process involves oscillations in the production of 
protein U, and that the concentration of U is given hj u = Uq sin out. In order to couple the 
oscillations of U to our network, we imagine inserting the gene encoding repressor adjacent 
to the gene encoding U. Then, since U is being transcribed periodically, the co-transcription 
of repressor will lead to an oscillating source term in Eq. (p!5D, 

i = m^fix) - - ixyxy + T sm{ujt) (16) 
y = rnyf{x)--fyy 

We first consider parameter values as in Fig. 2, so that the concentrations x and y oscillate 
in the absence of driving. Here, we are interested in how the drive affects the "internal" 
oscillations. Although there are many interesting properties associated with driven nonlinear 
equations such as Eq. (p!6D, we focus on the conditions whereby the periodic drive can cause 
the dynamics to shift the internal frequency and entrain to the external drive frequency uj. 
We utilize the numerical bifurcation and continuation package CONT |^9[ to determine the 
boundaries of the major resonance regions. These boundaries are depicted in the parameter- 
space plot of Fig. 5A, where the period of the drive is plotted versus the drive amplitude. 
The resonance regions form the so-called Arnold tongues, which display an increasing range 
of the locking period as the amplitude of drive is increased. Without the periodic drive, the 
period of the autonomous oscillations is equal to 14.6 minutes. As one might expect, the 
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dominant Arnold tongue is found around this autonomous period. Within this resonance 
region, the period of the oscillations is entrained, and is equal to the external periodic force. 
The second largest region of frequency locking occurs for periods of forcing which are close 
to half of the period of the autonomous oscillations. As a result of the periodic driving, 
we observe 1:2 locking, whereby the system responds with one oscillatory cycle, while the 
drive has undergone two cycles. Other depicted resonant regions (3:2, 2:1, 5:2, 3:1) display 
significantly narrower ranges for locking periods. This suggests that higher order frequency 
locking will be less common and probably unstable in the presence of noise. Outside the 
resonance regions shown in Fig. 5A one can find a rich structure of very narrow M:N locking 
regions with M and N quite large, together with quasiperiodic oscillations. The order of 



resonances along the drive period axis is given by the Farey sequence i.e. in between 
two resonance regions characterized by rational numbers, Mi:Ni and M2:N2, there is a region 
with ratio (Mi+M2):(Ni+N2). 

The preceding notions correspond to the driving of genetic networks which are intrinsi- 
cally oscillating. We now turn to a network designed with parameter values just outside the 
oscillatory region, and consider the use of resonance in the following application. Suppose 
there is a cellular process that depends critically on oscillations of a given amplitude. We 
seek a strategy for modifying the amplitude of this process if, for some reason, it is too small. 
For concreteness, consider a cellular process linked to the cell-division period of the host for 
our synthetic network. For E. coli cells at a temperature of ~ 37 degrees C, this period is 
of order 35 — 40 minutes. Using Figs. 3A and 4A, we can deduce parameter values that will 
cause a CI-RscA network to oscillate, when driven, with this period. The lower wedge of 
Fig. 3A implies that, for = 0.1, we should design the network with values of 7^. and jy 
just below the lower boundary of the wedge. Fig. 4A implies that, for = 0.1, a choice 
of 7y = 0.004 will yield oscillations with a period close to the cell-division period. In order 
to stay outside the oscillatory region, we therefore choose 'jy just below this value. Taken 
together, these choices will yield a network whereby oscillations can be induced by cellular 
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processes related to cell division. In Fig. 6A, we plot the drive versus response amplitudes (F 
vs F^) obtained from numerical integration of Eq. (|l^). We see that, depending on the prox- 
imity to the oscillatory region, oscillations are triggered when the drive reaches some critical 
amplitude. In Fig. 6B, we plot the gain g = (T + F^.)/F as a function of the drive amplitude 
F, and observe that, for certain values of the amplitude of the drive, the network can induce 
a significant gain. 



7 Harnessing the Lambda Switch 

The ability to switch between multiple stable states is a critical first step towards sophis- 
ticated cellular control schemes. Nonlinearities giving rise to two stable states suggest the 
possibility of using these states as digital signals to be processed in cellular-level computa- 



tions (see, for example, |3^, 0). One may eventually be able to produce systems in which 
sequences of such switching events are combined to control gene expression in complex ways. 
In any such application, the speed with which systems make transitions between their sta- 
ble states will act as a limiting factor on the time scales at which cellular events may be 
controlled. In this section, we describe a bistable switch based on the mechanism used by A 
phage, and show that such a system offers rapid switching times. 

The genetic network of A phage switches its host bacterium from the dormant lysogenous 



state to the lytic growth state in roughly twenty minutes 0]). As discussed in Section 2, the 
regulatory network implementing this exceptionally fast switch has two main features: two 
proteins (CI and Cro) compete directly for access to promoter sites; and one of the proteins 
(CI) positively regulates its own level of transcription. Here, we compare a synthetic switch 
based on the A phage's switching mechanism to another two-protein switch (the toggle switch 
described in Ref. [^), and numerically show that the A-like system offers a faster switching 
time under comparable conditions. 

To implement the synthetic A switch, we use the plasmid described in Section 3, on 
which the Prm promoter controls the expression of the A repressor protein, CI. To this, 
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we add a second plasmid on which the Pr promoter is used to control the expression of 
Cro. The operator regions ORl, 0R2, and 0R3 exist on each plasmid, and both proteins 
are capable of binding to these regions on either of the plasmids. On the PijAf-promoter 
plasmid, transcription of CI takes place whenever there is no protein (of either type) bound 
to 0R3; when CI is bound to 0R2, the rate of CI transcription is enhanced. On the Pr- 
promoter plasmid, Cro is transcribed only when operator site 0R3 is either clear, or has a 
Cro dimer bound to it; either protein being bound to either ORl or 0R2 has the effect of 
halting the transcription of Cro. 

Letting y represent the concentration of Cro, the competition for operator sites leads to 
equations of the form x = f{x,y) — 'y^x, y = g{x,y) — ■jyy. We derive the form of these 
equations by following the process described in Section 3. As with the CI plasmid of that 
section, we have Eq. 1 describing the equilibrium reactions for the binding of CI to the 
various operator sites. To these, we add the reactions entailing the binding of Cro, and the 
reactions in which both proteins are bound simultaneously to different operator sites: 



Y + Y 




Y2 


D + Y2 




Dl 


Dl + Y2 




DiD\ 


Dl + F2 


P2KA 


DlD\ 


DlDl + ^2 




DlD\D 


D^Df + Y2 


PaKa 


D\D^D 


D^ + Y2 


I35KA 





where Y represents the Cro monomer, and Df represents binding of protein p to the OKi 
site. For the operator region of A phage, we have /?i ~ /?2 — /^a ~ 0.08, and ~ /?5 ~ 1 

The transcriptional processes are as follows. Transcription of repressor takes place when 
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there is no protein (of either type) bound to 0R3. When repressor is bound to 0R2, the 
rate of repressor transcription is enhanced, and Cro is transcribed only when 0R3 is either 
vacant, or has a Cro dimer bound to it. If either repressor or Cro is bound to either OKI 
or 0R2, the production of Cro is halted. These processes, along with degradation, yield the 
following irreversible reactions. 



D + P ^ D + P + n^X (18) 

+ P ^ + P + n^X 

D^D^ + P "-^^ D^D^+P + n,X 

D + P ^ D + P + UyY 



Dl +P ^ Dl + P + UyY 



X 

Following the rate equation formulation of section 3, we obtain 



Q{x,y) 

where 



X = 7^7 ^ IxX 19 

Q{x,y) 
- - - lyy 



Q{x, y)^l + x' + aix^ + (JicrsX*^ + y' + (A + /Sa)^ + PiP3y^ + criPixV + P^xY- 
The derivatives are with respect to dimensionless time, with scaling as in Section 3; i 



t ( ktxPo drt^x V-^i -^2 ) , where ktx is the transcription rate constant for CI, and Ux is the number 
of CI monomers per mRNA transcript. The integers rrix and niy represent the plasmid copy 
numbers for the two species; Py is a constant related to the scaling of y relative to x. The 
parameters jx and jy are directly proportional to the decay rates of CI and Cro, respectively; 



24 



we will tune these values to cause transitions between stable states. The system exhibits 
bistability over a wide range of parameter values, and we plot the null-clines in Fig. 7A. 
For comparison, we now consider the co-repressive toggle switch briefly reviewed in sec- 



tion 2 This switch uses the CI and Lac proteins, where each protein shuts off transcrip- 
tion from the other protein's promoter region. The experimental design was guided by the 
model equations, 

ii = - u 20 









1 + 


u 




(l+l™)" 









where u and v are dimensionless concentrations of the Lac and CI proteins, respectively, 
and the time derivatives are with respect to a dimensionless time: r = /c^t, with = 
2.52h~^ p!2|, [T8| being the protein decay rate. The dimensionless parameters ai, 0^2, 6, and 



/i deflne the basic model. The CI protein used in the experiments is temperature-sensitive, 
changing its rate of degradation with temperature [p6|, |6^ ; we modify the original model 



slightly to include the factor k, which represents a varying decay rate for the CI protein. 
Switching is induced either by changing k, or by adjusting the concentration of isopropyl-/?- 
D-thiogalactopyranoside (IPTG); the parameters K = 2.9618 x 10~^ M and rj = 2.0015 deflne 
the effect of the inducer molecule IPTG on the Lac protein. Over a wide range of parameter 
values, the system has two stable flxed points; the null-clines are shown in Figure 7B. 

The time courses of switching between stable states in the two models are shown in 
Figure 8; transitions are induced by eliminating the bistability, then restoring it. The precise 
time course of switching from one stable state to another is determined by the way in which 
the model parameters are adjusted to eliminate the bistability. In each case, some parameter 
is increased until the system passes through a saddle node bifurcation: two stable flxed points 
and one unstable flxed point collapse into a single stable point. In an effort to examine the 
behaviour of the two systems under analogous conditions, we eliminate the bistability in 
every case by setting the system to a parameter value 10% past the bifurcation point. 
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The transitions shown in Figure 8 are generated as follows. The system begins (0-1 
hour) in its default bistable state, sitting at one of the two stable fixed points. Then (1-4 
hours) the bistability is eliminated (as described above), with the only remaining fixed point 
being such that the other protein has a high concentration. Once the concentrations have 
switched, the default parameters are restored and the system moves to the nearby stable 
fixed point (4-7 hours). Finally, the system is rendered monostable again (7-10 hours), 
causing another transition, followed by a period (10-11 hours) during which the bistable 
parameters are restored. 

Under the conditions shown, the A switch model displays significantly more rapid tran- 
sitions between its stable states than those seen in the model of the toggle switch. The 
numerical results indicate that the properties of the A switch do offer an advantage in terms 
of the speed of transitions, indicating that it may be fruitful to study synthetic models based 
on this natural system. Future analytical work on models such as the one presented in this 
section may allow us to make more precise statements regarding the source of this advantage. 



8 Conclusion 



From an engineering perspective, the control of cellular function through the design and 
manipulation of gene regulatory networks is an intriguing possibility. Current examples of 
potential applicability range from the use of genetically engineered microorganisms for envi- 
ronmental cleanup purposes [Q, to the flipping of genetic switches in mammalian neuronal 



cells While the experimental techniques employed in studies of this nature are certainly 
impressive, it is clear that reliable theoretical tools would be of enormous value. On a strictly 
practical level, such techniques could potentially reduce the degree of "trial-and-error" ex- 
perimentation. More importantly, computational and theoretical approaches will lead to 
testable predictions regarding the current understanding of complex biological networks. 
While other studies have centered on certain aspects of naturally-occurring genetic regu- 



latory networks P, 0, [1^, |2J, |25|, ^ an alternative approach is to focus on the 
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design of synthetic networks. Such an engineering-based approach has significant technolog- 
ical implications, and will lead, in a complementary fashion, to an enhanced understanding 
of biological design principles. In this work, we have shown how several synthetic networks 
can be designed from the genetic machinery of the virus A phage. We have highlighted some 
of the possible behavior of these networks through the discussion of the design of two types 
of switches and a relaxation oscillator. Additionally, in the case of the oscillator, we have 
coupled the network to an existing cellular process. Such coupling could lead to possible 
strategies for entraining or inducing network oscillations in cellular protein levels, and prove 
useful in the design of networks that interact with cellular processes that require precise 
timing. 

With regard to model formulation, there are several intriguing areas for further work. 
For one, the number of molecules governing the biochemistry of genetic networks is of- 
ten relatively small, leading to interesting issues involving internal noise. Recent pivotal 
work [l^, ^ has led to a systematic modeling approach which utilizes a Monte Carlo- 
type simulation of the biochemical reactions While this approach is impressively com- 
plete, its complexity makes analysis nearly impossible. An alternative approach could entail 
the use of Langevin equations, whereby the effects of internal noise are incorporated into 
stochastic terms whose magnitudes are concentration-dependent. Indeed, in the context of 
genetic switches, this approach has recently been suggested [Q. The advantage of this 
formulation is that stochastic effects can be viewed as a perturbation to the deterministic 
picture, so that analytic tools can be utilized. 

A potentially important technical issue involves the implicit assumption that the reactions 
take place in three-dimensional space. While this assumption is perhaps the most natural. 



proteins have been observed sliding along a DNA molecule in search of a promoter region [ 54 



so that protein-DNA reactions might effectively take place on a surface. While this would 
not alter the qualitative form of Eq. (BI), the exponents on the variable x could take on other 



values |65|, and this, in turn, could lead to significant quantitative differences. 
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It has been nearly 30 years since the pioneering theoretical work on interacting genetic 
networks Due, in part, to the inherent complexity of regulatory networks, the true 

significance of these studies had to await technological advances. Current progress in the 
study of both naturally occurring and synthetic genetic networks suggests that, as the pio- 
neers envisioned, tools from nonlinear dynamics and statistical physics will play important 
roles in the description and manipulation of the dynamics underlying cellular control. 
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Figure Captions 

FIG. 1. Results for additive noise with parameter value m = 1. (A) Bifurcation plot for 
the steady-state concentration of repressor versus the model parameter (B) The energy 
landscape. Stable equilibrium values of Eq. (|1^) (with D = 0) correspond to the valleys 
at z = —1.6 and 0.5, with an unstable value aX z = —0.52. (B) Steady-state probability 
distributions for noise strengths of = 0.04 (solid line) and D = 0.4 (dotted line). (C) 
The steady-state equilibrium value of z plotted versus noise strength. The corresponding 
concentration will increas as the noise causes the upper state of (B) to become increasingly 
populated. (D) Simulation of Eqs. (||) and (P) demonstrating the utilization of external 
noise for protein switching. Initially, the concentration begins at a level of [GFP] ~ 0.4 fiM 
corresponding to a low noise value of D = 0.01. After 40 minutes, a large 2-minute noise 
pulse of strength D = 1.0 is used to drive the concentration to ~ 2.2 fiM. Following this 
pulse, the noise is returned to its original value. At 80 minutes, a smaller 10-minute noise 
pulse of strength D = 0.1 is used to return the concentration to near its original value. The 



simulation technique is that of Ref. [63 . 



FIG. 2. The relaxation oscillator. (A) Schematic of the circuit. The Prm promoter is 
used on two plasmids to control the production of repressor (X) and RcsA (Y). After dimer- 
ization, repressor acts to turn on both plasmids through its interaction at Prm- As its 
promoter is activated, RcsA concentrations rise, leading to an induced reduction of repres- 
sor. (B) Simulation of Eqs. ([T5|) . Oscillations arise as the RcsA-induced degradation of 
repressor causes a transversal of the hysteresis diagram in Fig. lA. The parameter values are 
rrix = 10, my = 1, 7^ = 0.1, jy = 0.01, and j^y = 0.1. 

FIG. 3. Oscillatory regimes for the relaxation oscillator. (A) The bifurcation wedge is 
larger for smaller values of the parameter This larger regime corresponds to larger 

values of the RcsA degradation parameter 7^. Note that the native (i.e., without tuning) 
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degradation rates of 7^ ~ ~ 0.005 are very near the oscillatory regime. (B) Bifurcation 
diagrams as a function of and ^^^y^ and for two fixed values of ^y. The oscillatory regime 
is increased for smaller values of 7^, and, in both cases, small values of are preferable 
for oscillations. (C) The bifurcation diagram as a function of the copy numbers and 
rriy^ and for fixed degradation rates. Importantly, one can adjust the periodic regime to 
account for the unknown parameter 'y^^y. The figure also indicates that, for oscillations, one 
should choose as large a copy number as possible for the plasmid containing the repressor 
protein (ma,). In (A) and (B), constant parameter values are — 10 and my — 1, and in 
(C) 7a: = 1.0 and 7^, = 0.01. 

FIG. 4. Parameter dependence of the oscillatory period. (A) An increase in 7^ decreases 
the period of oscillations. (B) The period depends very weakly on the copy number. In (A) 
— 10 and in (B) — 0.01, and for both plots, other parameter values are 7^, = 0.1, 
IxV — 0.1, and my — 1. 

FIG. 5. Dynamics of a periodically driven relaxation oscillator - Eq.(15); — 0.1, 7^ = 0.01, 
Ixy — 0.1. (A) Resonant regions in period-amplitude parametric plane. Solid lines (limit 
lines of periodic solutions) together with dashed (period doubling) lines define boundaries 
of stable periodic solutions for a given phase locking region M:N, where M is the number of 
relaxation oscillation and N is the number of driving sinusoidal oscillations. (B,C,D) Oscilla- 
tions in periodically driven repressor (top curve) concentration together with the oscillations 
of the sinusoidal driving (bottom curve). (B) 1:1 synchronization; the 14.6 minute period of 
cl oscillations is equal to the driving period. (C) 1:2 phase locking; the 29.2 minute period 
of the cl oscillations is twice as long as the driving period. (D) 2:1 phase locking; the 7.6 
minute period of the cl oscillations is equal to one half of the driving period. 

FIG. 6. (A) As a function of the driving amplitude F, the amplitude F^; of the induced 
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network oscillations shows a sharp increase for a critical value of the drive. The critical 
value corresponds to a drive large enough to induce the hysteretic oscillations, and it in- 
creases as one decreases 7j, and moves away from the oscillatory region in parameter space. 
The three curves denoted 1,2, and 3 are for •jy values of 0.0038, 0.0036, and 0.0034. (B) The 
gain as a function of the drive amplitude for 7^, = 0.0038. Close to the oscillatory region, a 
significant gain in the drive amplitude can be induced. Parameter values for both plots are 
7a; = 0.1, 'jxy = 0.1, rrix = 10, and niy = 1. Note that, corresponding to these values, the 
network does not oscillate (without driving) for 7^ < 0.004 (see the bottom wedge of Fig. 3A). 

FIG. 7 Null-clines for the two-protein bistable switch systems. Stable fixed points are marked 
with circles, and unstable fixed points are marked with squares. (A) Null-clines for the syn- 
thetic A switch, Eqs. (|T9|) . Solid line: a; = cline. Dashed line: y = cline. Parameter 
values: = 0.004; -fy = 0.008; Py = 62.92; a = 11; = rriy = 1; (Xi = 2; (T2 = 0.08; 

= = /53 = 0.08; and (3^ = (j^ = 1. (B) Null-clines for the toggle switch, Eqs. (ID. Solid 
line: ii = cline. Dashed line: v = cline. Parameter values (from Ref. p2|): ai = 156.25; 
02 = 15.6; S = 2.5; /i = 1; = 2.0015; [IPTG] = 0; k = 1. 

FIG. 8. Transitions between stable states for the two-protein bistable switch systems. The 
protein concentrations have been normalized (the trace for each protein is normalized rela- 
tive to its own maximum value). The system parameters are varied over time, altering the 
stability of the system and causing transitions, as described in the text. Upper Plots: The 
switching of Lac (solid) and CI (dashed) in the synthetic toggle model ||2^ . The parameter 
values are as given in the caption to Fig. 7, except as follows. (1-4 hours): [IPTG] = 2 
mM, k = 1.0. (7-10 hours): [IPTG] = 0.0, k = 50.81. Lower Plots: The switching of 
CRO (dashed line) and CI (solid) in the synthetic A model. The parameter values are as 
given in the caption to Fig. 7, except as follows. (1-4 hours): 'jx = 0.004, ■jy = 21.6. (7-10 
hours): 7^. = 18.0, 7^ = 0.008. 
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Figure 1 - Hasty et al. 
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Figure 3 - Hasty et al. 
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Figure 4 - Hasty et al. 




Figure 5 - Hasty et al. 




0.0 0.2 0.4 0.6 0.8 1.0 

G 



g 



30 



20 



10 








0.0 0.2 0.4 0.6 0.8 1.0 



Figure 6 - Hasty et al. 




Figure 7 - Hasty et al. 




Figure 8 - Hasty et al. 



